DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS 


TECHNICAL  NOTE  2748 

ON  TRANSONIC  FLOW  PAST  A  WAVE -SHAPED  WALL 

By  Carl  Kaplan 

Langley  Aeronautical  Laboratory 
Langley  Field,  Va. 


/ 


Washington 
August  1952 


20000508  214 


IT 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  27^8 

ON  TRANSONIC  FLOW  PAST  A  WAVE-SHAPED  WALL 
By  Carl  Kaplan 


SUMMARY 


The  present  paper  is  an  extension  of  a  previous  investigation 
(described  in  NACA  TN  2383)  concerned  with  the  solution  of  the  nonlinear 
differential  equation  for  transonic  flow  past  a  wavy  wall.  In  the 
present  work  several  new  notions  are  introduced  which  permit  the  solu¬ 
tion  of  the  recursion  formulas  arising  from  the  method  of  integration 
in  series.  In  addition,  a  novel  numerical  test  of  convergence,  applied 
to  the  power  series  (in  the  transonic  similarity  parameter)  representing 
the  local  Mach  number  distribution  at  the  boundary,  indicates  that  smooth 
symmetrical  potential  flow  past  the  wavy  wall  is  no  longer  possible  once 
the  critical  value  of  the  stream  Mach  number  has  been  exceeded. 


INTRODUCTION 


In  NACA  TN  2383  (ref.  1)  the  writer  considered  the  problem  of  two- 
dimensional  transonic  flow  past  an  infinitely  long  sinusoidal  solid 
boundary.  The  problem  was  treated  in  the  physical  plane  and  the  purpose 
was  to  investigate  whether  or  not  the  flow  past  the  wavy  wall  remains  a 
smooth  symmetrical  type  of  potential  flow  when  the  undisturbed-stream 
Mach  number  exceeds  its  critical  value.  By  a  smooth  type  of  potential 
flow  is  meant  one  for  which  the  velocity  potential,  say,  and  its  first 
derivatives  are  single-valued  and  continuous;  that  is,  there  are  no 
discontinuities  of  the  nature  of  shock  waves. 


L 


[nitially,  the  Prandtl-Busemann  small-perturbation  method  was 
applied  and  the  velocity  potential  developed  inclusive  of  the  third 

power  in  the  "thickness"  coefficient  6  =  — 


A 


where 

a  amplitude  of  wavy  wall 

wave  length  of  wavy  wall 


X 


2 
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The  velocity  potential  was  then  referred  to  the  critical  speed  of  sound 

k(l  -  M2)3//g 
7  +  1 


ccr,  the  coefficient  e  was  replaced  by 


where 


k  = 


(7  +  l)e 


is  the  transonic  similarity  parameter,  and  terms 


(1  -  M2)3/2 

involving  powers  of  1  -  higher  than  the  first  were  neglected 


(main  assumption  for  transonic  flow) .  This  simplified  or  transonic 
form  of  the  Prandtl-Busemann  solution  was  shown  to  be  identical  with 
the  one  obtained  by  solving  the  simplified  nonlinear  differential  equa¬ 
tion  for  transonic  flow,  with  the  boundary  condition  taken  not  at  the 
wave-shaped  wall  but  at  the  flat  plate  corresponding  to  vanishing 
amplitude-  The  calculation  was  carried  through  the  sixth  power  in  the 
transonic  similarity  parameter  k  and  corresponds  to  the  insuperable 
task  of  obtaining  the  Prandtl-Busemann  solution  to  the  sixth  power  in 
the  thickness  coefficient  e-  Thus  each  iteration  step  of  the  Prandtl- 
Busemann  method  contributes  to  the  transonic  form  of  the  solution,  which 
may  therefore  be  considered  a  result  of  thin-profile  theory  with  disturb 
ances  not  necessarily  small  compared  with  1  -  Moo^.  The  main  conclusion 
reached  in  reference  1  was  that  the  transonic  similarity  parameter  k 
must  be  less  than  4/3  —  a  value  still  somewhat  greater  than  the  critical 
value  0.83770  calculated  there. 


The  purpose  of  the  present  work  is  to  express  the  solution  of  the 
problem  of  transonic  flow  past  the  wavy  wall  in  a  form  more  suitable  for 
general  considerations  and  to  prove  that  the  assumed  smooth  symmetrical 
type  of  potential  flow  cannot  exist  at  stream  Mach  numbers  beyond  the 
critical  value. 


SYMBOIB 

x,  y  nondimens ional  rectangular  Cartesian  coordinates 

a  amplitude  of  wavy  wall  \ 

X  wave  length  of  wavy  wall 

X/2«  reference  element  of  length 

£L 

e  "thickness"  coefficient  of  wavy  wall, 

7  ratio  of  specific  heats  at  constant  pressure  and  constant 

volume 

U  stream  velocity 
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M  local  Mach  number 

Moo  undisturbed-stream  Mach  number 

( 

k  transonic  similarity  parameter. 


ccr 

critical  speed  of  sound 

0 

velocity  potential  of  flow 

f(x,y) 

transonic  disturbance  potential 

fn(y) 

functions  of  y  only,  related  to  f(x,y) 

An  P 
q  t 

numerical  coefficients 

An  P 
q 

generating  functions  of  k, 

00 

E 

r=0 

An  P  lc2r 

Aq  p+r  K 

•^m  p 

functions  of  dummy  variable 

r. 

00 

V  n  p  n 

L ^n-m  r  >  l°wer  label 

V  ii  iu  ' 

n=l,m 

starts  from  1  when  m  is  negative  and  from  m  when  m  is 
positive 


Primes  denote  differentiation  with  respect  to  independent  variable. 

ANALYSIS 

General  Formulas 

If  the  undisturbed  stream  is  in  the  direction  of  the  positive  x-axis, 
then  the  velocity  potential  0  referred  to  the  critical  velocity  ccr 
can  be  written  as  (see  ref.  1) 


(7  +  l)e 

“^372 


(l  -  M2)' 


0  =  x  +  y  +  2(1  ~  Moo2)f(x,y) 
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where  the  second  term  on  the  right-hand  side  is  a  disturbance  velocity- 
potential  and  implies  that  terms  involving  powers  of  1  -  Mro2  higher 
than  the  first  have  been  neglected.  The  differential  equation  for 
transonic  flow  satisfied  by  the  function  f(x,y)  is  obtained  from  the 
general  differential  equation  for  compressible  flow  and  takes  the 
following  simplified  nonlinear  form: 


d2f  _  df  dff 

dy2  Sx  8x2 


The  boundary  conditions  to  be  fulfilled  by  f(x^y)  are  as  follows: 


> 


(at  y  =  °°) 


r  (2) 


=  -k  sin  x  .  ■  (at  y  =  0,  <  x  <  ») 

dy 


Here  x  and  y  are  nondimens ional  rectangular  Cartesian  coordinates 
simply  related  to  the  physical  plane  coordinates  X  and  Y  by  means  of 
the  trans format ion 


f- 


hi. 


x  =  X 

y  =  (l  -  M*2)1/2! 

and  the  equation  of  the  wavy  wall  is  Y  =  €  cos  X  or 

y  .=  (i  _  M002)1//2e  cos  x.  Clearly,  f(x,y)  involves  only  the  variables 
x  and  y  and  the  transonic  similarity  parameter  k. 

The  most  general  form  for  the  function  f(x,y)  to  ensure  symmetrical 
potential  flow  past  the  wavy  wall  is  the  following  (see  fig.  l): 
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f(x,y) 


-x 


+ 


00 

r 


fn  sin  nx 


(3) 


where  the  fn  are  functions  of  y  only.  When  this  form  for  f(x,y) 
is  substituted  into  equation  (l)  and  the  coefficients  of  the  separate 
harmonic  terms  sin  nx  are  equated  to  zero,  the  following  system  of 
second-order  nonlinear  ordinary  differential  equations  for  fn  results: 


n-1 


Ln 


n2fn  = 


1 

¥ 


nLm(n-  m)fmfn-m  “  |  n  H  m(n  +  m)fmf. 


m=l 


m=l 


n+m 


(n  =  1,  2,  .  .  .  oo)  (4) 

Before  proceeding  to  the  solution  of  these  equations,  several  formulas 
of  general  interest  and  subsequent  use  are  given.  They  have  been  derived 
in  reference  1.  The  local  Mach  number  distribution  is  given  by 

1  -  M2  =  -(l  -  Mro2)|£  (5) 


H 


The  equation  from  which  the  critical  value  of  the  transonic  simi¬ 
larity  parameter  is  calculated  follows  from  equation  (5)  by  taking 
x  =  0,  y  =  0,  and  M  =  1;  that  is, 


=  0 


or  with  the  aid  of  equation  (3), 


L_  nfn(0)  =  1  (6) 

n=l 


A 


6 


'v,K 


p  - 


_22 


The  pressure  coefficient 
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is  given  by 


_2_( 

7  +  1' 


')(>  *  g) 


Integration  in  Series  for  the  Functions  .  fn 

In  reference  1,  equations  (4)  were  solved  by  an  iteration  pro¬ 
cedure  and  approximate  expressions  for  fp  to  fg  were  obtained.  An 

examination  of  these  expressions  showed  that  the  general  form  of  fn  is 

n-1 

2p+n-2  oo_ 

tn=Ie-(2^  H 

p=0  q=0  r=p 

(n  =  1,  2,  .  .  .  °°)  (8) 


where,  if  p  =  0,  the  upper  limit  of  q  is  n-1  and,  if  P  /  0,  the 
upper  limit  of  q  is  2p  +  n  -  2.  The  four-labeled  coefficients  A£  f 

are  real  numbers  calculated  from  recursion  formulas  obtained  from  the 
system  of  differential  equations  (4)  and  the  boundary  condition  at  the 
surface  of  the  wavy  wall.  The  boundary  condition  at  y  =  <»  is  auto¬ 
matically  satisfied  by  the  form  of  fnj  whereas  the  boundary  condition 
at  the  wall  takes  the  form 


^  - -  -  S  /  *  y.  __ 

The  four-labeled  coefficients  A£  P 


(fn’)y=0  - 


(n  =  1) 
(n  {  1) 


Inserting  the  expression  for  fn  given  by  equation  (8)  into  equations  (9) 
yields  immediately  the  following  results: 


(2P  +  n)A£ 

P=0 


p=0,l 


An  P 
A1  r 


(n  =  1,  2, 
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where,  if  n  =  1,  the  lower  limit  of  p  on  the  right-hand  side  is  unity 
and,  if  n  ^  1,  the  lower  limit  of  p  is  zero.  Also,  if  n  =  1,  the 
upper  limit  r  of  p  goes  from  1  to  «  and,  if  n  f  1,  r  goes 
from  0  to  oo. 


By  elementary  manipulation  of  series,  the  second  of  equations  (10) 
can  he  replaced  by  the  following  more  useful  forms: 


nAn  0  =  An  0 


0  0 


1  0 


(n  =  2,  3*  •  •  •  °°) 


0 

r 


p=0,l  P=1 


(r  =  1,  2, 

Vn  =  1#  2, 


(ID 


where  in  the  first  term  on  the  right-hand  side  the  lower  label  p  starts 
from  1  when  n  =  1  and  from  0  when  n  /  1. 


In  reference  1,  recursion  formulas  were  derived  for  the  coef¬ 
ficients  A^  p.  In  the  present  paper  a  much  more  significant  approach 

is  introduced.  Note  that  equation  (8)  can  be  rewritten  in  the  following 
form: 


00 

1"  / 


(12) 


where 


00 


An  P 
q  p+r 


k2r 


In  a  manner  similar  to  that  described  in  reference  1  for  the  coef¬ 
ficients  A^  p,  recursion  formulas  can  be  obtained  for  the  power 

series  An^P.  Indeed,  the  two  types  of  recursion  formulas  are  intimately 

connected  and,  for  a  given  value  of  p,  the  one  can  be  obtained  from  the 
other  by  mere  inspection.  A  single  recursion  formula  can  be  written  for 
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the  general  quantity  Anqp  but  the  resulting  expression  is  cumbersome 

and  serves  no  practical  purpose.  It  is  much  more  desirable  to  obtain 
separate  recursion  formulas  for  each  value  of  p.  As  examples  of  pro¬ 
cedure,  the  recursion  formulas  for  p  =  0,  1,  and  2  are  considered  in 
the  sections  which  follow. 


Recursion  Formula  for<  A 

With  p  =  0,  the  recursion  formula  is  (compare  with  eq.  (57)  of 
ref.  1) 


/  .  n  0 

2n(q  +  l)Aq+1 


.  .  ,  .  n-2  n  0 

(q  +  l) ( q  +  2)5  q  Aq+2  + 


\ -  (m  +  l)(n  -  m  -  DA™'1  0  Am+1  0 

Z _  L _  q-i  <3-^1 

qx=0  m=q-q1 

(n  =  2,  3>  •  •  .  00)  q  =  0,  1,  .  .  .  n  -  2)  (13) 


where 


n-2  0 


(q  =  n  -  2) 

(q  ^  n  -  2) 


This  recursion  formula  can  be  solved,  the  solution  starting  with 
q  =  n  -  2  and  descending  towards  q  =  0.  Thus,  for  q  =  n  -  2,  equa¬ 
tion  (13)  becomes 


n/  ,  N„n  0  'v —  t  ,  w  n\,m+l  0  .n-m-1  0 

8(n  "  l)An-l  =  }_  (m  +  l)(n  "  111  '  1)A  m  A  n-m-2 

m=0 

(n  =  2,  3,  •  •  •  °°)  (14) 

Wow,  multiply  both  sides  of  this  equation  by  rn  where  r  is  a  dummy 

variable  and  sum  from  n  =  2  to  n  =  oo.  Then 


oo  oo  n-^ 

Y~  (n  -  1)aJJ  i  rn  =  V~  rn  Y~  (m  +  l)(n  -  m  -  l)A] 

*  1  7  /->  _ 


m+1  0  ^n-m-1  0 
m  n-m-2 


n=2  m=0 
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then 


\  _  V~  An  0 

l10  -  Z_  n-l  r 


A1  o' 


=  y  nAn  ®  rn-l 
L —  n-l 


°°  n-l 

(A1  O')2  =  )  r11"2  F  m(n  -  m)Am  0  An"m  0 

i— -  i—  m-1  n-m-1 


By  observation,  it  can  be  seen  that  the  right-hand  side  of  equation  (15) 
is  equal  to  (rAx  Q'j2.  Equation  (15)  can  thus  be  replaced  by  the 
following  first-order  nonlinear  differential  equation: 


(rAi  o')2  -  8rAl  0’  +  8Ai  0  =  0  (■ 

The  solution  of  this  equation  is 

irA  * 

A  »  _  r  JTa1  0 

•A1  0  “  coe 

where  cQ  is  the  arbitrary  constant  of  integration.  From  the  defini¬ 
tion  of  Aj_  0  it  follows  that,,  with  r  =  0, 


c  =  A1.0 
co  R  0 


and  hence 


A1  o' 


^rA1  0  Al  °' 

0  ^  0 


(17) 


10 


NACA  TN  27bQ 


Note  that  by  definition  An  ®  involves  the  coefficients  A*?  ®  which 

are  ultimately  calculated  by  means  of  the  auxiliary  relations,  equa¬ 
tions  ( 11) ,  engendered  by  the  boundary  condition  at  the  surface  of  the 
wavy  wall.  In  solutions  of  the  recursion  formulas,  therefore,  the 

coefficients  A11^  appear  as  undetermined  quantities. 


A 


The  expansion  of 


1  0 
,1  0 


in  powers  of  rA' 


1  0 


is  a  nontrivial 


problem  which  fortunately  can  be  solved  with  the  aid  of  Lagrange's 
investigations  on  the  reversion  of  power  series.  In  reference  2, 
Lagrange's  problem  is  illustrated  by  the  following  example: 


To  expand  eax  in  powers  of  y 


xe 


bx 


The  result  is  expressed  as 


,ax 


1  +  ay  + 


i(a  -  2b)  2  ,  a(a  -  31) 2  1 


2! 


3! 


+ 


+ 


a(a  -  rib) 11-1  n  + 


(18) 


In  particular,  with  a  =  1,  b  =  -1,  and  |  =  e,  the  solution  of  the 
equation  log  |  =  y|  or  |  =  e^  (generally  referred  to  as  Eisenstein's 
problem)  is  given  by  the  series 


n11-2 


L-  (n  -  i): 
n=l 


7  y 


n-l 


Then,  if  |  is  replaced  by 
of  equation  (17)  is 


12  and  y  by  ^  rA^^,  the  solution 


A 


1  0 
0 


A  1 

A1  0 


n-2 


n 


A' 


1  0 


n=l 


(n  -  D'U11-1^  ° 


aA  °f 1  n-i 
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or 


nAn  0  r11"1 

Z_  n-1 


00 


n 


n-2 


(n  -  l) !  l+n-1 


n 

r 


n-1 


Hence,  by  equating  coefficients  of  equal  powers  of  r  on  both  sides  of 
this  equation, 


n  0 
'  -1 


(n  =  2,  3)  •  •  •  “) 


(19) 


which  is  the  solution  of  the  recursion  formula,  equation  (l4).  The 
coefficient  A^-^  can  be  considered  as  the  generating  function  for  the 
set  of  coefficients  A^^  Thus,  equation  (19)  can  be  written  as 


OO 


A 


n 

n-1 


0  k2r 
r  ■K 


1  0 
0  r 


(20) 


Then,  by  equating  coefficients  of  equal  powers  of  k  on  both  sides  of 

this  equation  and  putting  A^  ®  =  1  (see  boundary  condition,  eq.  (10)), 
the  following  equations  are  obtained: 


n  0  _  nn-2 
An-1  0 


n!  4 


n-1 


,  n  0 
^-1  1 


A  n  0 

^n-l  2 


nn-2 


a; 


1  0 


(n-l)!!*"-1  01 


n 


n-2 


Ai  °  + 


n 


n-2 


(n  -  1)  ! 4n_1  0  2  2 !  ( n  -  2) ! 4: 


s)‘ 


(21) 


12 
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2  0  1 

Note  that  from  the  first  of  these  formulas,  A  =  —  and  that  from 

10  o 

the  boundary  conditions,  equations  (ll),  Aq  q  = 

Consider  now  q  =  n  -  3*  Then  equation  ( 13)  becomes 


IX  1 

2n(n  -  2)A^_°  =  (n  -  l)(n.  -  2)A^_J  +  |  n  m(n  -  ^n-m-1  ^22) 

m=2 

f 

» 

Multiply  both  sides  of  this  equation  by  rn  and  sum  from  n  =  3  to 
n  =  00.  Then 


2  L  n(n  -  2)a£_°  rn  =  V  (n  -  l)(n  -  2)A^  r11  + 

n=2  n=l 


n  0 


1 1: £ -u - a- °  (23) 


n=3  m=2 


Let 


*1 0  -  r  cl 


n=l 


and 


A2  0  =  H  An-2  rn 
n=2 


then 


00  n 


ll  o'A2  0*  =  H  rI1  n  m(n  -  m  +  2)^-1  An“n- 


n=l  m=l 


m  0  ,n-m+2  0 
m 


NACA  TN  2rjbQ 


13 


t 


U 


It  can  be  shown  easily  that  the  last  term  on  the  right-hand  side  of 

equation  (23)  is  equal  to  —  r^r^A^  o*^2  0*)'*  •®erice^  equation  (23) 
can  be  rewritten  as  the  following  differential  equation: 


rA 


1  0 


'(rAl  o')'  =  o" 


Then  equation  (24)  can  be  written  as 


(24) 


2r%  Q"  -  2rA2  0'  =  r%  0”  -  2rAx  0'  +  2AX  0  +  I  r(r\  Q'A2  y)' 

Wow,  from  equation  (l 6),  it  follows  that 

(rAi  o')2  =  8(rAl  O'  -  A1  o) 

(rAi  o') '  ‘  ~-AL°  p.  '>  <®> 


Let  rAg  q’  = 
Then  equation 


(2 6)  becomes 


where 


v  is  the  new  dependent  variable. 


-  6rAl  0 


*  +  4a 
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with  the  aid  of  equations  (25), 


■■  •  ih 


Hence 


T  '  e  rAi  o'  +  5  ci 


8rA 2  o' 


(rAi  o'  )3  ,  „  (rAi  o') 

rnw +  ci 


The  arhritary  constant  c-^  is  determined  by  differentiating  this  last 
equation  twice  and  evaluating  for  r  =  0.  Thus, 


Finally, 


=  6h  - *-r- 

(AV)' 


A2  0'  =Kr2|l0’)5  +  8  ^  °" 


Wow,  from  equation  (l8)  with  a  =  2,  b  =  -1,  and  I  -  ex,  it  follows 
that 


nn~3  n-2 
(n  -  2)  ! 
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Or  with 


A  ' 

A1  0 


> •  i  «v 


(■>  ”')'■■  isrrS^M-'" 


Mi  ■  ’  £ 


Thus,  equation  (27)  becomes 


Y  <  °  r"-1 .  |  f-  — p(a1o°)>-1 
A-  n‘2  2  Z_  (n  -  3):4n‘2V  °  ' 


,2  0  “ 
A  0  V 


I>“  - 1)A-? rn_1 


Hence,  replacing  A^_^  by  its  value  from  equation  (19)  and  equating 

the  coefficients  of  equal  powers  of  r  on  both  sides  of  this  equation 
yields 


An  o  =  8  _ P*1.'.3  . 

n-2  °  .  , 

(n  -  2)lk 


3  (a1  0)n' 

„,n-lVA  0  ) 


2  0  +  1  n11-4  /  1  0\n 

°  2(n-  3)!4n-A  °> 


(n  -  3*  •  •  .oo)  (28) 
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the  solution  of  the  recursion  formula,  equation  (22).  Again,  equa¬ 
tion  (28)  can  be  written  as 


00 

J2  An-2  ?  k 


noo  .00  \  n-2 

(n-2)!4nltr0  [ho  / 


n11'4  Ar~  a1  0  k2r\n 

tTn1!  /  0  r 

’  3>:4  \feo  ) 


Then,  equating  the  coefficients  of  equal  powers  of  k  on  both  sides  of 
this  equation  and  putting  aJ®  =  1  and  -^o  0  =  lS^  yields1 


ln-2  0 


2(n-2)!kn“1  2(n-3)!)+n"1 


.  n  0  8nn“3  ,2  0  ,  n11"3  fll  0 

ln-2  1  "  “  ,  n-1  A01  ,  .  n-1  0  1 

n  -  2) ! 4  (n  -  3)54 


.  n  0  _  8nn~3  ,20  nn~3 

n-2  2  ,  .  .  n-1  0  2  Q  % , >  • 

(n  -  2) 14  (n  -  3) 


_n-3 

— T  A 


A1  °  +  8A1  0  A2  0  + 
A0  2  0  1  0  1 


l/.l  0\21  .  n11"3  f,l  OV 

l) J  +  2(n  .  4>!^A°  ^ 


3  0  1 

Rote  that  from  the  first  of  these  formulas  q  =  2% 

3  0  1 

the  boundary  conditions,  equations  (ll),  Aq  ©  =  yry 


and  that  from 


3T 
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The  procedure  followed  in  order  to  obtain  equation  (28)  for  A^ 
is  a  general  one  and  with  very  little  difficulty  other  members  of  the 
family  An  0  can  be  obtained. 

For  example. 


From  this  equation  it  follows  that 


n  0  _  2  nn-4 


n-3  0 


1  n1' 


3(n-3)!4n~1  2  (n  -  4)!4n~1 


a  n  0  |,0  n11-^  ,3  0  0  n11-^  ,2  0 

An-3 1  -  ^  ; - — nn  A0 1 +  8 ; — — nn  Ao  1 + 

(n  -  3)  !4  (n  -  4)  !4 


nn~k  .1  0[3  1  1  1  ~| 

L.n-1  0  1_3  (n  -  4)!  2  (n  -  5) 


From  the  first  of  these  formulas  and  the  boundary  conditions,  equa- 
tions  (11),  aJ  °  =  3^  and  a£  g  =  5^5. 


At  this  point,  it  is  noted  that  the  coefficients  of  the  form  A^  q 

are  calculated  from  the  first  formula  of  each  set,  equations  (21),  (30), 
(32),  and  so  forth.  A  number  of  this  type  of  coefficient  have  been 
evaluated  (see  ref.  l)  by  means  of  the  recursion  formulas  for  the 


l8 
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coefficients  A^  £  themselves.  They  are  listed  as  follows: 


.2  0 

1 

A6  0  _ 

7  x  13 

A1  0 

■  8 

A10' 

72  x  256" 

.3  0 

1 

a7  0  _ 

13 

Ai  0 

=  25 

Ai  0  " 

18  x  256 

A  0 

7 

A8  0  - 

13  X  19 

A1  0 

38^ 

A1  0 

576  X  258 

50 _ 

■1  0  "  768 


A  careful  examination  of  these  numerical  values  leads  to  the  general 
rule. 


where  by  definition 

|3n  _  5}  =  1  x  k  X  7  X  10  X  13  X  .  .  .  X  (3n  -  5) 

In  the  expression  for  the  local  Mach  number  distribution  evaluated 
at  the  crest  of  the  wavy  wall  (x  =  0,  y  =  0),  there  occurs  the  following 
power  series: 


00 


NACA  TN  2748 


19 


F  =  k  + 


L,  n!4n_1 


This  power  series  can  he  expressed  in  the  closed  form 


/  3  \2/3 

F  =  2  -  2^1  -  kj 

The  graph  of  F  against  k  is  a  semicuhical  parabola  with  the  cusp 

4 

point  at  k  =  —  and  F  =  2.  With  the  necessary  condition  that  one  and 

only  one  value  of  k  correspond  to  a  given  value  of  F,  the  transonic 
similarity  parameter  k  cannot  be  greater  than  4/3.  Moreover,  the 
lower  limit  of  k  is  zero  when  the  amplitude  of  the  wavy  wall  is  zero 
but  Moo  Is  different  from  unity. 


Recursion  Formula  for  A 


With  p  =  1,  the  recursion  formula  for  A11^  is  (compare  with 
eq.  (69)  of  ref.  l) 

4(n  +  l)Anq1  -  2(n  +  2)(q  +  1)5^  +  (q  +  l)(q  +  2)5n“^n  A^  = 

"  I  n6^  H  L  m(n  "  ffl)An"^i°  Aq-qi  "  2  n(n  +  1)Alo°  An+q  ° 

q1=0  m=q-q1 


(n  =  1,  2,  .  .  .  00)  q  =  0,  1, 


n)  (33) 


where 


5n-l,n  _  0 

q  1 


( q  =  n  -  1  or  n) 
(q  i  n  -  1  or  n) 


_  0 
6q  -  1 


(q  =  n) 
(q  /  n) 
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The  solution  of  this  recursion  formula  proceeds  as  in  the  case  p  -  , 
the  solution  starting  with  q  =  n  and  descending  towards  q  =  0.  Thus, 
for  q  =  n,  equation  (33)  becomes 


1  .10  An+1  0 

ft  nA  0  A  n 


(n  =  1,  2,  •  •  •  00 ) 


CV  -  -  ¥*  -  1)AV  An-“ 


(n  =  2,  3, 


Hence,  inserting  the  expression  for  aJJ_J  given  by  equation  (19)  gives 

An-1  1  =  _  i - IL-1 - -(a1  °)n  (n  =  2,  3,  •  •  •  “)  (3^) 

n-1  8  /  _  0\ ,  I,11-1  \  0  / 


(n  -  2)'X 


f-.n-l  1  ,  Pr  _  1  n»-3  f  r  A1  0  V2r\n+1 

11-1  r+1  8  (n  -  °  r  / 

From  this  relation,  by  equating  coefficients  of  equal  powers  of  k  on 
both  sides  of  the  equation,  the  following  equations  are  obtained 

An-1  1  .  .  i - - - 

n-i  1  8  (n  _  2) ! hn_1 


,  n-1  1 


n-1  2  o  (n  _  2) 


.n-1  1 
n-1  3 


3  n11’3  ,10.1  ^3  A1  0 

o'  :  rmr  Ao  1  «  ,  ,..i.n-i  o  1 


14n_1  01  8  (n  .  3)!^ 


3 

n-3  r 

1  0  +  (a1  0>1 

21  -  i  — 

nn-3 

=  ‘  8  (n 

-  zy.k^L 

•0  2  +  ^A0  1/ 

J  8  (n 

-  3)!4nTIl 

„  / . 1  o\ 

2]  1 

n11"3  /nl 

°\2 

3(A0  ij 

-  4)!^n_1^  ° 

1/ 

(n  -  2,  3,  •  •  •  “) 
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Consider  now  q  =  n  -  1)  equation  (33)  becomes 


4(n+l)A^  =  2n(n+2)Ann1  -  \  n  Y1  (m  +  l)  (n  -  m  -  l)Am+ffl1  0  A11^"^1  - 


,1  0  An+1  0 


JL  /  ,  -  \  .  JL  U  A  li'TX  V. 

-  n(n  +  l)A  Q  A  n_1 


(n  =  1,  2,  .  .  .  °o) 


Multiply  both  sides  of  this  equation  by  rn  and  sum  from  n  =  1  to 
n  =  oo.  Then 


oo  oo 

4  (n  +  rn  =  2  n(n  +  rn  - 


i  r  nrn  £  (m  +  1)  (n  -  m  -  l)A"+±  0  a""-1,1 


n=l  m=0 


m  n-m-1 


OO 

1  A1  0  YL  n(n  +  l)An+1  ’ 

2  0  n  n-1 


0  rn 


00  „  -1 

A  =  A  1  r11 
A0  1  Z_  A  n  r 
n=l 


A  _  V~  fln  1  -r11 
A11'L  n-1  r 

n=l 


i  o  -  r  c; rn 

n=l 


-  r  A"  |  r“ 


2  0  / —  n-2 

n=2 


Then  it  can  easily  be  shown  that  the  second  term  on  the  right-hand  side 
of  equation  (36)  is  -  i  r^r^A^  i'Ai  o')'  an<^  "^at  equation  (36)  can 
be  replaced  by 


22 
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4(rAi  i)'  =  2T\  i  +  6rAo  i*  "  \  r(r2Ao  l'Al  0*)' 


-  rA1  0  A  " 
2  0  2  0 


(37) 


Now, 


^  1  *  -  g  A®  Y  <n  -  X>A!L?  rn_1 

n=l 


or,  with  the  aid  of  equations  (25), 


rAor‘^  a1o°K  o’) 

Hence,  from  the  use  of  equations  (25)  and  equation  (27)  for  A2  q'  it 
follows,  after  routine  calculations,  that  the  integral  of  equation  (37) 
is 

rAi  1  =  lis  a1o°(-  l(rAi  o’) 2  “  s(rAi  o’)3  "  i(rAi  o’)4  + 


(ci  +  h)  |(rAi  o')2  "  2r2Ai  0. 


3} 


where 


.  2  0 


c^  =  64 


K°y 


Now  from  equation  (l8), 

( 


m 


rAn  '  \  =  m 

1  0 


f  S™±-(Ay)n  rn 


(38) 
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Then 


Equating  coefficients  of  equal  powers  of  r  on  both  sides  of  this 
equation  gives 


(n  =  2,  3,  .  .  .  °°)  (4-0) 

From  this  equation  the  expressions  for  all  the  coefficients  of  the 
n  1  1 

type  r+-L  with  n  =  2,  3>  •  •  •  00  and  r  =  0,  1,  .  .  .  oo  can 

be  obtained.  Thus,  for  example, 


2k 
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n-1  1  _  1  nn-5r_20  ,  70 

An-2  1  32  411-I [(n  -2)'.  (n  -  3)! 


37 

+  (n-_  Vr  +  i 


n  5)0 


,n-l  1  _  _1_  n 


n-5 


n-2  2  32  411-I 


(5n  +  1)aJ  J  +  3£A§  fW^aJT  + 


(5n  +  1)Aq  °  +  32Aq  ?](„  ^3yf  +  jpTu  +  5)aJ  “  + 
®6Ao  D(n  -\-)T  8A0  l)(„  -5):] 


>  (41) 


For  q  =  n  -  2,  equation  (33)  becomes 


Mn  +  l)An  J  =  2(n  -  l)(n  +  2)A^  -  n(n  -  ljA11^- 


n  1 


>n  1 


n-2 


n-2 


1  \  ,  W  ^ .m+2  0  . n-m-2  1 

i  n  >  (m  +  2)(n  -  m  -  2)A  m  A  n_m_2  + 

d  mR)  L 

,  '  w  ,  um+l  0  .n-m-1  1| 

(m  +  l)(n  -  m  -  1)A  m  A  n_m_2  J  ~ 


1  (  .  tuI  0  An+1  0 

2  n<n  +  1)A  0  A  n-2 


(^2) 


The  solution  of  this  recursion  formula  follows  along  the  same  lines  as 
that  for  equation  (36)  and  leads  to  the  following  result: 


4T 
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From  this  equation  there  follow,  in  the  usual  manner,  expressions  for 
coefficients  of  the  type  A^  Thus,  for  example,  with  r  =  0, 


n-3  1 


/8l  ,  21  2  1  3\  1  (\ 

k  16  n  +Zn)w-~3T'.'+  V 


201  +  3  n2  + 

6k  2  n  + 


'205  9_ 

535  32 


+ 


+ 


j£L  1  nn~7 
256  (n  -  6)!  i^n-1 


m 
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Recursion  Formula  for  A 


With  p  =  2,  the  recursion  formula  for  An  2  is 


8(n  +  2)An  2  -  2(n  +  4)(q  +  l)8n+2  An  2  +  (q  +  l)(q  +  2)&n+1'ri+2  A^J  = 
q  q  q+1  y  ^ 


n+i-qi 

1  n6n+2  sa-l,n,n+l  80,1,2  p  V-  (m  -  2)  ( 

2  q  1 1  m  ^ — 

1  q^O  m=q-q1 


n  -  m  + 


n-m+2  0  m-2  2  _  1  ^n+l,n+2  n  X  1+  m(n  .  1  Am  1 


q  =0  m=q-q1 


1  ,  t  _,/Rn+2  .1  0  .n+1  1  .  Rn+l,n+2  .1  1  fin+l  0 

-n(n  +  1)^5  q  A  0  A  g  +  S  '  A  0  A  q  + 

0  n+2  1  1  n+1  o\  /  „\ /cn+2  .2  0  An+2  0  .0  .2  0  .n+2  o\ 

8°'q  Al^l  ^+1^  .  n(n  +  2)  ^8  q  A  Q  A  q  +  8q  A  ±  A  ^  j 

(n  =  1,  2,  .  .  .  oo>  q  =  0,  1,  .  .  .n  +  2)  (45) 

where  delta  is  defined  in  the  usual  manner.  The  solution  of  this 
recursion  formula  proceeds  as  in  the  previous  cases,  starting  with 
q  =  n  +  2  and  descending  towards  q  =  0.  Thus,  for  q  =  n  +  2,  equa¬ 
tion  (45)  becomes: 


-  i  nA2,0  Am!,» 
R  1  n+1 


(n  =  1,  2,  .  .  .  oo) 


,n-2  2 


jkn  -  2)A21°  A^_°  (n  =  3>  4,  .  .  .  °°) 
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Then,  from  equation  (±9) , 


or 


(n  =  3,  4,  . 


2 

r+2 


k2r 


0$  (46) 


By  equating  coefficients  of  equal  powers  of  k  on  both  sides  of  this 
equation,  the  following  equations  are  obtained: 


An'2 

n 

i 

ii 

CM  OJ 

(n  -  2)nn“2 

•  i.n+2 
n!4 

An-2 

2 

(n2  -  4)nn“2 

n 

3 

,  .,n+2 
n:4 

An-2 

2 

(n2  -  4)nn_2 

n 

4  "  “ 

i  i,n+2 
n!4 

0 

1 


> 


(47) 


Consider  now  q  =  n  +  1$  equation  (45)  becomes 

8(n  +  2)A^+2  =  2(n  +  2)(n  +  4)A^+2  - 

_n-2 

1  ^  ,  ,  , w  , s .m+1  0  .n-m-1  2 

2nL(m+  l)(n  -  m  -  DA  A  n_m+1 
c  m=0 

i  n(n  +  1)(a10°  A”^1  +  A1^  An+J  °)  - 
n(n  +  2)(a2q°  A^0  +  A2-^  An+2  °) 


( n.  —  lj  2}  •  •  •  oo)  (48) 
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With  the  introduction  of 


uu 

-  Z  Cp  ^ 


-2  2  fes.  tl+2 


00 

V"  n  < 

Ag  o  -  l-2  An-2  rl 


„  v-  An  2  n 

A_1  2  =  2_  Vl  r 

n=l 


A  =  y  An  1  r11 
K0  1  l—  A  n 
n=l 


lio  =  t 


An  ?  r11 


equation  (48)  is  replaced  by  the  following  ordinary  linear  differential 
equation: 

8(r2A_x  2) '  =  2t2(tA_2  2')  '  +  12r(rA_2  2'j  +  l6rA_2  2  - 

|  r2(r2A!  0«A_2  2‘)  -  §  A1Q°  r2AQ  ±”  -  l  A1,1  r2A,  - 


2  1  10 


a20°  rAi  o"  +  A  0  A1  o'  -  A  1  rA2  o"  + 


A2  0  A  .  _  A1  0  A2  ° 
120  0*0 


By  repeated  use  of  equations  (25 )>  the  solution  of  this  differential 
equation  is  found  to  he  as  follows: 


8r2A-l  2  =  (a10°)2  i(rAl  o’)  -  i(rAl  o')2  -  768 (rAl  o' 


i^(rAl  o')4  +  k  r2Al  o"  +  \  A2  o  +  ^(cl  +  k)  H(rAl  of 
o]}  -  r(A10°)3  +  A2o°|^1  0  -  |(rAlVoj^  - 


rA1  0  A2  0 
0  0 


(50) 
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Finally,  with  the  aid  of  equation  (38)  and  the  definition  of  A_-j_ 


From  this  equation,  there  follows  in  the  usual  manner  formulas  for  the 
coefficients  of  the  type  aJ“^  .  Thus,  for  r  =  0, 


DISCUSSION  OF  CONVERGENCE  OF  SMOOTH  TYPE  OF  POTENTIAL 
FLOW  PAST  WAVY  WALT. 


In  the  preceding  sections,  a  number  of  examples  of  recursion 
formulas  and  their  solutions  have  been  given  in  considerable  detail. 

The  purpose  of  this  exposition  is  threefold:  First,  to  show  the  inherent 
elegance  of  the  method  of  integration  in  series  although  the  equations 
concerned  are  nonlinear  in  character)  second,  to  present  a  type  of 
analysis  which  may  be  useful  in  other  problems  involving  nonlinear 
differential  equations;  and  third,  to  indicate  that  an  analytical  proof 
of  convergence  may  ultimately  be  obtained  by  careful  examination  of  the 

recursion  formulas  for  the  quantities  An  ^  and  their  solutions.  One 

M 


example  is  the  obtaining  of  the  general  expression  for  the  coefficients 
n  0  j, 

A1  0  and  the  subsequent  conclusion  that  k  ^ 


In  actual  practice  it  has  been  found  more  convenient  to  evaluate 
the  coefficients  A^  ^  from  their  separate  recursion  formulas  rather 
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"than  to  derive  the  general  formulas .  The  appendix  contains  the  exact 
numerical  values  of  the  coefficients  necessary  for  the  development  of 
the  functions  fn  to  the  eighth  power  in  the  transonic  similarity 

parameter  k.  These  coefficients  are  utilized  to  demonstrate  numerically 
the  test  for  convergence  of  the  smooth  symmetrical  type  of  potential  flow 
(eq.  (3))  assumed  in  this  paper.  Thus,  when  the  form  of  fn  given  by 
equation  (8)  is  inserted  into  equation  (5)  for  the  local  Mach  number 
distribution  and  evaluated  at  the  surface  of  the  sinusoidal  wall  (y  =  0), 
the  following  result  is  obtained: 


M2  -  1 

2 

1  -  M*, 


-1  + 


kn  (n  “  2m)  cos(n  - 
n=l  m=0 


m 

2m)  x  ^ 

p=0 


An-2m 

0 


P 

m 


(52) 


where  ^  denotes  the  integral  part  of  n/2.  At  the  crest  of  the  wavy 

wall  (x  =  0) ,  the  point  of  maximum  fluid  velocity,  with  the  numerical 
values  for  the  coefficients  inserted,  equation  (52)  becomes: 


M2  -  1 

1  >  V 


359381  .6  .  7396757  k7  +  81688733  k8  + 

’  '  g  K  2  2 

270  x  256  648o  x  256  864oo  x  256 


The  critical  value  of  k  (that  is,  when  M  =  l)  calculated  from  this 
equation  is 

kcr  =  0.83244 


(Note  that  in  ref.  1  the  value  kcr  =  0.83770  corresponds  to  the  first 
six  terms  of  eq.  (53).)  Consider  now  the  infinite  series 


where 


An0,  JlLli} 

0  0  nn!4n-1 


(n  =  2,  3,  •  •  •  °°) 
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An  0 

The  Cauchy  ratio  test  Rln  =  ^  yields  in  the  limit  n->«  the 

A  0  0 

result  that  the  radius  of  convergence  Rj  is  equal  to  or  less  than  I4-/3 . 

If  the  corresponding  ratios  R2n  are  formed  for  the  right-hand  side  of 

R 

equation  (53)  and  the  quotient  Nn  =  —  is  calculated,  the  resulting 

■^Pn 

sequence  of  numbers  is  as  follows: 


n 

Rln 

R2n 

Rn 

1 

16.0 

8.0 

2.0 

2 

1.92 

2.34375 

3 

3.04762 

1 . 78041 

1.7H75 

4 

2.5 

1.33366 

1.87472 

5 

2.21546 

1.35000 

l.64io4 

6 

2.04167 

1.17721 

1.73432 

7 

1.92481 

1.19588 

1.60950 

The  noteworthy  feature  of  this  table'  is  that  although 
sumably  Rpn)  Is  converging  quite  slowly  toward  Rp  = 


Rln  Pre_ 

^  (and  R2),  the 


quotient  Wn  exhibits  a  strong  tendency  to  approach  an  asymptotic 

value  N  for  a  relatively  low  value  of  n.  Figure  2  shows  this  tendency 
in  a  graphic  manner.  The  apparent  asymptote  represented  by  the  straight 
line  is  the  ratio  of  4/3,  the  limit  of  Rln  as  n— *00,  and  of  0.83244, 


the  critical  value  of  k.  Certainly,  the  rapid  approach  of  the  lower 
dotted  curve  toward  the  apparent  asymptote  and  the  decreasing  oscilla¬ 
tions  represented  by  the  upper  dotted  curve  indicate  that  the  critical 
value  of  k  is  the  radius  of  convergence  R2  of  the  power  series  on 
the  right-hand  side  of  equation  (53).  Thus,  the  critical  value  of  the 
stream  Mach  number  marks  the  limit  of  convergence  of  the  smooth  sym¬ 
metrical  type  of  potential  flow  assumed.  The  ability  to  approximate 
closely  the  limiting  value  N  is  a  matter  of  luck)  namely,  the  choice 
of  the  known  comparison  series.  Once,  however,  the  proper  comparison 
series  has  been  selected  and  the  approach  to  an  asymptote  indicated,  there 
can  be  no  question  of  the  meaningfulness  of  the  approximate  value  of  N. 
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obtained.  It  may  be  that  one  would  like  to  extend  figure  2  to  n  =  9- 
This  extension  would  entail  the  forbidding  calculation  of  an  additional 

185  coefficients  An  P.  The  .result  presumably  would  be  to  decrease 

q  r 

slightly  the  approximate  critical  value  of  k  and  thereby  raise  slightly 
the  straight-line  asumptote  of  figure  2.  This  extension  of  figure  2 
would  show  still  more  convincingly  the  approach  to  an  asymptote  and  the 
dying-out  of  the  oscillations.  Perhaps  more  important  still,  figure  2 
definitely  shows  that  conclusions  based  on  less  than  six  or  eight  terms 
are  mere  speculations  in  this  field. 


CONCLUDING  REMARKS 


If  the  numerical  test  of  convergence  presented  is  acceptable,  the 
conclusion  to  be  drawn  is  that  smooth  symmetrical  potential  flow  past 
the  wavy  wall  exists  only  for  the  purely  subsonic  range.  Moreover  no 
such  flow  can  possibly  represent  the  transonic  or  mixed  type  for  which 
a  local  region  of  supersonic  flow  near  the  solid  boundary  is  imbedded 
in  the  otherwise  subsonic  stream.  It  follows  as  a  corollary  that  the 
transonic  or  mixed  type  of  flow  past  the  wavy  wall  is  necessarily  an 
asymmetric  one.  This  asymmetry  in  the  flow  pattern  entails  a  resistance 
usually  defined  as  wave  drag.  As  shown  by  experimental  observations,  the 
shock  wave  associated  with  wave  drag  closes  the  downstream  portion  of 
the  local  supersonic  zone. 

As  a  final  remark  —  although  the  analysis  and  conclusions  of  the 
present  work  refer  directly  to  the  wavy  wall,  the  suggested  result  that 
the  critical  stream  Mach  number  marks  the  limit  of  smooth  potential  flow 
very  likely  applies  to  other  boundaries.  This  conclusion  is  based  on 
the  idea  that  the  gradual  transition  from  a  purely  sinusoidal  wall  to  a 
boundary  composed  of  a  single  bump,  say,  introduces  no  essential  changes 
in  the  analysis. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  May  6,  1952. 
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APPENDIX 


NUMERICAL  VALUES  OF  THE  COEFFICIENTS  An  P 

q  r 


Coefficients  for  f^: 


A 


1  0 
0  0 


=  1 


A1  1  =  .  _5_ 

oi  256 


A1  2 
ft0  2 


_ 65 _ 

576  x  256 


.1  3  _  3365 

0  o  0 

1152  X  256 


A1  0  =  JL 

0  1  256 


A1  1  =  -  2765  A1  2  =  26435  A1  3  =  3907 

02  9  x  2562  03  81  X  2563  1  3  432  x  256s 


1  0  =  1861  ,11..  184345  ,12..  7 

0  2  3  x  2562  3  1728  x  256s  1  2  72x256 


Al  3  =  £3 

23  12  x  2562 


A10=  4896755 
03  81  x  2563 


A1  1 
A1 1 


1 

54 


a1  2 
1  3 


53995 

1728  x  256s 


a1  3  =  .  583 

33  54  x  2562 


A 


1  1  =  33 

12  '  64  x  256 


A1  2 
2  2 


9 

32  x  256 


Af  3  .  .  _ 119 _ 

^  3  12  X  2562 


A1  1  =  .  139 

1  3  4  x  2562 


A1  2  =  .  £765 

23  96  x  2562 


A1  3  - 
?  3 


1 

96  x  256 


12  1 

13  2  '  6  X  256 


A1  2  =  _  55 

33  8  x  2562 
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Coefficients  for  fg'- 


A 


2  0  _  1 

oo  15 


A0  1 


125 _ 
36  X  256 


.2  2  =  12245 

02  144  X  256s 


.2  3  _  .  17999 

03  90  x  2563 


A 


2  0  _  419 

0  1  72  x  256 


.2  1  =  _  45895 

0  2  108  x  2562 


a2  2  =  39023983 
0  3  1620  x  2563 


a2  3  =  72977 

1  3  5760  x  2562 


20=  234215 

02  432  x  2562 


2  1  _  4210409 

03  81  x  256s 


a2  2  _  6245 

12  72  x  2562 


a2  3  =  45409 

2  3  2160  x  2562 


a2  0  =  42791533 
03  648  x  256s 


k2  1 
ftl  1 


a2  2  =  104063 
13  90  x  2563 


a2  3  =  3919 . 

33  864  x  256s 


A 


2  0 
1  0 


1 

8 


a2  1  =  .  815  _ 

12  576  x  256 


A2  2  =  _  _ 79 _ 

2  2  3  x  296s 


A2  3  _  _  ,...iag22 _ 

^  3  864  x  2562 


A 


2  0 
1  1 


_  ll. 
1024 


2  1  =  456185  A2  2  =  _  167 

1  3  3456  x  256s  23  4  x  256s 


A2  3  =  .  73 _ 

53  8  x  2562 


A2  0  _  4085 

1  2  24  X  256s 


A2  1  =  -  JL 
2  1  128 


.2  2  _  47 

32  "  192  x  256 


25 

12  X  2562 


2  0  _  21287  A2  1  =  -  11 

1  3  324  X  2562  22  32  X  256 


ft2  2 

k3  3 


241 99 

864  x  256s 


4  3 


4811 


192  x  256 


A2  2 
r4  2 


.2  2 
&4  3 


12  x  256 


11 


2  x  256 
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35 


a3  0  =  139 

22  2  x  2562 


36 
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Coefficients  for  f^: 

- 

.4  0  _ 

7 

11 

1 — 1 1 — 1 

-4-  0 

c 

35251 

A4  2 

741115 

Ao  0  ' 

6  x  256 

90  X  2562 

A0  2 

.  5184  x  256s 

4- 

H  O 

II 

390547 

a4i  = 
0  2 

80872738 

A4  2 

38381086 

2 

720  x  256 

675  x  2563 

A1  2 

405  x  2563 

II 

O  OJ 

"=<p 

1704729613 

A4  1  = 
1  1 

2465 

A4  2 

423833 

16200  x  256 

576  x  256 

2  2 

1152  x  2562 

.4  0 

7 

11 

H  OJ 

-4-  H 
< 

30596538 

A4  2 
3  2 

20845 

A1  0  = 

384 

135  X  2563 

144  x  2562 

.4  0  _ 

32947 

Ak  1  = 
2  1 

455 

A4  2 

173 

Ai  1  " 

36  x  256s 

96  x  256 

a4  2 

12  x  256s 

20527213 

Ao  0  " 
2  2 

127825 

A4  2 

_  30 

A1  2  - 

162  x  2563 

192  x  2562 

5  2 

2562 

A4  °  _ 
2  0 

3 

AU  1  = 
3  1 

119 

A4  2 

■  6  2 

36 

128 

48  x  256 

5  x  256^ 

A1*'  0  = 
2  1 

617 

A4  1  - 

51343 

288  x  256 

A3  2  " 

216  x  256s 

a4  0  _ 
2  2" 

369587 

A4  1  - 

25 

1728  x  2562 

A4  1 

48  x  256 

a4  °  _ 

l 

A4  1  - 

275 

3  0 

96 

A4  2  " 

8  x  2562 

0  =  _11 _ 

3  1  24  X  256 


4  0  _  4811 

3  2  i44  x  256s 
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Coefficients  for  fc 


A5  0  =  7 

0  0  384o 


1 

0  1 


A5  0  _  90354659 

0  1  259200  x  2562 


1  1 


a5  0  =  7 

10  7^8 


A5  1 
A2  1 


a5  o 
A1  1 


886243 


1080  x  256' 


l5  l  = 
3  l 


'5  0  _  25 


2  0  1536 


>5  1. 
l4r 


A- 


5  0  _  24375 


2  1 


32  x  256" 


a5  1 
5  l 


\P  0  = 

3  0 


9 

384 


£  0  =  8950 

31  27  X  256s 


A5  0 

a4  0 


l5  0  = 

l4i 


25 


24  x  256 

1379 

24  X  256£ 


5470394 
81  X  256s 

-  88739 

96  x  2562 

25  x  15587 
288  x  2562 

25  X  997 

24  x  2562 

105 

64  x  256 

9 

32  x  256 
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Coefficients  for  fg: 


A, 


0 

91 

A6 

1 

2483629 

0  : 

432  X  258 

Ao 

1 

54  x  2563 

0  . 

285522407 

A6 

1  _ 

4371541 

1  ■ 

486o  x  2583 

1 

1 

5760  x  256s 

0  . 

91 

a6 

1  _ 

265901 

0  ' 

’  72  x  256 

2 

1 

192  x  2562 

0  - 

144862403 

A6 

1 

403535 

1 

810  x  2563 

3 

l 

288  x  2562 

0 

23 

A6 

l  _ 

79841 

0  " 

'  8  x  256 

l 

96  x  256^ 

0 

2553091 

A6 

l  _ 

2205 

1  ' 

2880  x  2562 

A5 

1  “ 

8  x  2562 

0  . 

.  13 

A6 

1  _ 

2401 

0  ' 

1024 

6 

l 

60  x  2562 

A6  0  =  ^387 

3  1  96  x  2562 


0 

0 


15 

8  x  256 


a6  0  =  819 

^  1  4  x  256s 

a6  0  -  9 

5  0  20  x  256 

A6  0  =  297 

51  10  x  2562 
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Coefficients  for  fy: 


a7  0  =  _ 13_ . 

00  126  X  256 


A7  0  =  13 

10  18  X  256 


A7  0  _  287 

2  0  i44  x  256 


A7  0  =  833 

30  288  X  256 


A7  0  _  343 

4  0  i44  x  256 


A7  0  _  343 

50  320  x  256 


A7  0  =  2401 

60  2 

45  X  256 


4o 
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Coefficients  for  fg: 


.8  0  _  13  X  19 
00  18  x  2562 


a8 .0  =  13  x  19 
10  576  x  256 


.8  0  _  89 

20  64  x  256 


A8  0  =  59 

30  24  X  256 


a8  0  „  ■lt2_ 

4  0  18  x  256 


A8  0  =  — i2— 
5  0  2880 


A8  0  =  _L_ 
6  0  2880 


A8  0  -  _ i — 

7  0  2520 


NACA  TN  2rjk8 


hi 


REFERENCES 


1.  Kaplan,  Carl:  On  a  Solution  of  the  Nonlinear  Differential  Equation 

for  Transonic  Flow  Past  a  Wave-Shaped  Wall.  NACA  Rep.  1069,  1952. 
(Supersedes  NACA  TN  2383.) 

2.  Bromwich,  T.  J.  I‘a. :  An  Introduction  to  the  Theory  of  Infinite 

Series.  Second  ed.,  Macmillan  and  Co.,  Ltd.,  19^2,  pp.  l60-l6l. 


I 


CO 

TJ  7J  t- 

sc.h 

lua 

E  i2  "& 


Q  g> 

•  w  .2 

is1 

lit 


uj 

S**  2 
•2  <)  •<.. 

®  H  “ 

§P 

•gft< 

. 

533 

>»P«i  2< 
o  o  s2  S' 

f  5  B  7J  S 

pg  13  O  rt  N 


◄  igdg 

^ZO^£ 


)  »  a 
>  £  S'  o  >> 

1  p  a  **  a 

^  S  T3  U 

,  a>  S?  a)  efl 

1  **-<  71  S 

©  o  g  B 

^  tr  **h 


1 3  |  .2  « 


a  -a  5 

cn  > 

^  Q) 
cd  ^ 

c  £  0 


■  y  u  « 

M  C  B  ® 

!  5  "I  £>  S 


I  Z  >  Si 

:  s  §  s 


i  a  3.2 

;  §  g * 


C  66  .  ° 
On )  © 
03-5—1 

o  to  > 
e  rt  *-j 

8  a  3 

>  •'H 

£  a  £  B 
O  73  a>  -a 

Sf  a  3 

13  S  a  o 

5  4J  fl)  n 

©  o  S  ® 

B  « 

s!  I  & 

~  E  *a  .q 
Jh  3  3 

"  S  o  i 
®  a  a,  a 


00 

■a  7;  t- 

8  rt  « 

is  o  z 
SC.H 

te-J  < 

E  t2  J5 


Q  SP 

•  w  .2 

IS! 13 

2g 

A)  k  • 

c<S 
6-^2 
**-i  c 

®  H  3 

Sg  21 

g2c 

. 

u  SI 

►.h  & 

55  oO^  c? 

|  9  5  -s  S 

«  5  o  S  « 

!Z-<guZ 

Sl*dS 

Lc  IS  jz;  <  $ 

zzo&S. 


JC  V. 

12.2< 
j  8  ®  a 
5  O  H  ^ 

:s| b  3 

i  t*  S  S  ’S, 
I"  2|  8 
!•!  ft®  A 
1 2  ■a  S  « 

!  §  0)  «  o 
_  a  h  u, 


l-gl; 

rt  ” 

rt  ajs  , 


o  >>  o 

*M  -*J  © 

Tj  '5  a 

aa  “ 
a  a « 


-  5  ©  *-< 
PS  “35 
s  ■a  <  £  a 

3  b  «  .S  rt 

a  2  'S  ’3  P 

2  Z  «  a  f 

3  rt  “  £  ^ 

)  2  ^  (D 
*  2  —  xj  H 

«■§  5  ts 

3  ^  3  ■g  a 
!  a  h  S  53 
}.  g  of  "j  "2 

i  e  -s  1  5 
1  §  Z,  8  S 

'  rt  Qi  2 
1  £  o  £  a 


4)  fi  ^ 
1  §  8  § 
;  S  s  o 
2  S 
f  S  *  8 
1  §  J  ® 

i  8  7s 

!  °-S  » 

[“art 
1  •£  a>  o> 

^  S3  T3 

“«  5 

h  t,  " 
Oni  « 

a  I  & 

3  O  -S 
e  Pa  O 


0) 

I  2  m 

f  ® 

I1  a 

CQ  <*h 

&  o 

li 

> 

ED 

rt  73 

a  g 
^  ^3 
o  7;  . 

a  S-o 

73  ©  TJ 
”  £  0) 
■tt  53  © 

S  p  8 

■*7  © 

O  JZ  © 

a  >  co 

73  ©  t 

Z T  r— (  W 

•2  B  II 
'£  ®  a 

ft)  QQ  Q 

a  ai 


